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Abstract 

I study the Coulomb drag between two identical layers in the Anderson in- 
sulting state. The dependence of the trans-resistance pt on the localization 
length ^ is considered, for both a Mott insulator and an Efros-Shklovskii (ES) 
insulator. In the former, pt is monotonically increasing with in the latter, 
the presence of a Coulomb gap leads to an opposite result: pt is enhanced 
with a decreasing ^, with the same exponential factor as the single layer re- 
sistivity. This distinction reflects the relatively pronounced role of excited 
density fluctuations in the ES state, implied by the enhancement in the rate 
of hopping processes at low frequencies. The magnitude of drag is estimated 
for typical experimental parameters in the different cases. It is concluded 
that a measurement of drag can be used to distinguish between interacting 
and non-interacting insulating state. 
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I. INTRODUCTION AND PRINCIPAL RESULTS 



A rich variety of phenomena in disordered electronic systems is associated with the An- 
derson locahzation The fundamental effect is a manifestation of quantum coherence: a 
subtle destructive interference of multiple partial waves in a random medium leads to ex- 
ponentially localized electron eigenstates. In three-dimensional (3D) conductors, a metal- 
insulator transition (MIT) can be driven by tunning either the disorder potential or the 
carrier density across a threshold value. The transition is identified as second order: it is 
associated with a divergence of a correlation length ^, which in the insulating side has the 
physical significance of a localization length. A true phase-transition of this kind is generi- 
cally absent in two-dimensions (2D), at vanishing magnetic fields 0: instead, a 2D system 
exhibits a smooth cross-over between a "weak insulator" (in which at finite temperatures 
T the conductivity a acquires a negative logarithmic correction), and a "strong insulator" 
(characterized by an exponentially small cr(T)). The presence of a magnetic field as well 
as spin-orbit scattering ( 0, 0), dramatically alter this behavior - they suppress local- 
ization, and in principle can recover the metallic state in 2D. Most prominently, in strong 
magnetic fields, singular extended states play an important role in the mechanism for the 
quantum Hall (QH) effect [Q. Recent experimental studies have posed new challenges in 
the understanding of Anderson insulators. For example, in QH systems, the observation of 
a weak-to-strong localization cross-over near filling fraction 1/2, is possibly an interesting 
indication for B = Q localization of "composite Fermions" In contrast, a true MIT has 
been observed at i? = 0, whose origin is yet obscure. In all the above, electron-electron 
interactions and their interplay with the disorder |^ play an important role, and complicates 
the theoretical analysis. 

In view of the ongoing research activities described above, it is highly desired to have an 
extended array of different probing techniques. In the present paper, I suggest that Coulomb 
drag in a double layer system is potentially an interesting probe, that can diagnose subtle 
differences between distinct insulating states. 
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Coulomb drag is a manifestation of the coupling between two spatially separated 
systems of charge carriers, due to Coulomb interactions across a barrier separating them. 
In the presence of a transport current in one layer, density fluctuations in that layer exert 
a frictional force on the other, and consequently induce a voltage in the latter (in an open 
circuit configuration) - even when tunneling between the layers is suppressed. The strength 
of this effect is characterized by the trans-resistivity, pt = -E2/J1, where E2 is the parallel 
electric field induced in layer 2 in response to a current density ji established in layer 
1. Experimentally, drag has been observed in various semiconductor heterostructures [§; 
theoretically, it has been a subject of much recent activity |[T0|-[16|. In particular, it was 
suggested as a useful test of certain electronic states in the layers, such as compressible 
states in QH systems [l^Jl^ , and superfluid electron-hole condensates . 



As pointed out in Ref . |jT2[ , inter-layer drag can probe the dynamics and response prop- 
erties of electronic systems in a domain that is inaccessible to transport measurements in 
a single layer. As it stems from a frictional force. Coulomb drag is dominated by the in- 
teraction between relatively long-lived density fluctuations in the two layers. Therefore, its 
low temperature behavior is sensitive to the ability of the electronic systems to create and 
maintain such fluctuations, reflected by the density-density response functions xi,2{(l,^) at 
low frequencies {u — 0) and finite wave-vectors [|l^ q. To see this, note that to leading 



order in the screened inter-layer interaction U, the trans-resistivity can be expressed as |11 



P^"^ f _fq_ 2 p , I U (g, ^) pimxi (g, ^)Imx2 {q, uj) , . 

vm(i)n(2)e2i (27r)2Mo 4smh\f3hu;/2) ' ^ ^ 

here P = 1/kBT where T is the temperature, n'-^^n'-^^ are the carrier densities in layers 1,2, 
and 

V{q)e-'''^ 

^ [1 + V{q)xi{q, CO)] [1 + V{q)x2{q, uj)] - V^{q)e-^<^'^Xi{q, uj)X2{q, 00) 
where V{q) = 27re^/eg (e being the background dielectric constant), and d the inter-layer 
distance. In Eq. (|T|), the integration over u is effectively cutoff by T, while the cutoff on q 
is set by 1/d; the former is stricter at low T. 



In the vicinity of a MIT, the behavior of drag is expected to be particularly interesting, 
due to the significant role played by density fluctuations. Variations in the conductivity 
(of one layer or both) have competing effects on the availability of density fluctuations 
and on their stability. In the conducting state (e.g., in the diffusive regime [|11|]), the drag 
is enhanced as the diffusion coefficient is reduced, since then the decay rate of density 
fluctuations is slowed down. On the other hand, in the insulating state the creation of 
density fluctuations is suppressed at low q and u, and this can lead to a reduction of pt with 
a reduced localization length - provided that the simultaneous variations in the screening 
properties are not dominant |]l4| , p!8| . Thus, in contrast with the single layer resistivity, pt 
is potentially a non-monotonous function of the parameter which drives the transition; a 
similar behavior is expected also near a smooth (nonOOcritical) cross-over between a weak 
and strong insulator. At the transition, density fluctuations and their dynamics become 
critical, and one expects a pronounced peak in the drag (similarly to the behavior predicted 
in [|12| for QH transitions). 

In this paper I study the Coulomb drag in Anderson insulators, in order to verify under 
what circumstences the above described qualitative picture holds. I calculate pt in a double- 
layer system, where the electronic states in the layers are assumed for simplicity to be 
identical and uncorrelated. Two possible insulating states are considered: a Mott insulator 
||19|| , where the in-layer long-range Coulomb interactions are suppressed, and an Efros- 
Shklovskii (ES) state where the interactions are signiflcant. In both cases, the dominant 



transport mechanism is assumed to be variable range hopping ||21[] among localized sites, 
and Pt is evaluated as a function of the localization length ^. 

The resulting trans-resistance indicates a dramatic difference between the Mott and ES 
insulators. While the former exhibits a suppression of pt with a decreasing C,, as implied 
by the naive argument pointed out earlier, in the latter the dominant contribution to pt 
is enhanced with the same exponential factor as the single layer resistivity. As will be 
shown below, this follows from the relative enhancement of hopping processes at — 0, 
due to the presence of a Coulomb gap in the ES state. Drag measurement is therefore 



a suggestive experimental means of distinguishing the two types of insulating states. In 
contrast, the single layer resistivity exhibits a qualitatively similar dependence on ^ and T: 
p{T) ~ e'^'^o/T)°' ^ where in the 2D Mott insulator a = 1/3 and To ~ and in the ES state 
a = 1/2 and To ~ 

In the following sections I detail the calculation of pf. in Sec. fT|, for the Mott insulator, 



and in Sec. |T| for the ES insulator. The conclusions and experimental implications are 



summarized in Sec. IV. 



II. MOTT INSULATORS 

I consider two parallel 2D layers separated by a perfect insulating barrier of thickness 
d, in both of which the electronic state is a non-interacting Mott insulator. The trans- 
resistance is evaluated using Eqs. (|l]), @, where cu) = X2{(1,^) is the density response 
function, which has been derived diagramatically by Vollhardt and Wolfle ||22|. To leading 
order in q and u, it is given by 

dn Dq^ 

= dHW^^^^^TVi^^^^] ■ 

Here dn/dfi is the density of states, r is the ellastic mean free time, and D = VpT/2 [vp 
being the Fermi velocity) is the diffusion coefficient in the conducting state. The localization 
effect is represented by the frequency c^o: as pointed out in Ref. ||22|, Eq. is consistent 
with a simple hydrodynamical model, in which localization is manifested as an effective 
restoring force acting on density fluctuations, with a characteristic oscillator frequency uq. 
The localization length is directly related to luq through ^ = vp/\^iUQ. 

1 next assume that T, the effective upper cutoff on u in Eq. (|l|), is sufficiently low that 
the u <^ u!o,uJqT. The real and imaginary parts of xiQ,^) (hereon, the layers indices are 
dropped) are then given, respectively, by 

dn (7^ 

Imx(g,..)-|^-^,^^^^. (5) 
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Substitution in Eqs. (|T|) and (0) yields the final expression for pt. The result depends 
crucially on whether ^ is smaller or greater then the layers separation and exhibits a 
cross-over between the two limiting cases considered below. 

Deep in the insulating state, ^ <^ d, one can neglect the term in the denominators 
of Eqs. and (in the effective range of wave-vectors, q < 1/d and hence q -C 1/0- 
The inter-layer interaction can be then approximated by its unscreened form, U {q, u) ^ 
V"(g)e~'^'^. The integrations in Eq. (||) are streight-forward and result in 

where n = ra^^^ = n^'^\ and qrp = 2'Ke^{dn/ djj) / e is the Thomas-Fermi screening wave- 
vector. 

When the localization length is increased (e.g. by controlling a parameter that drives 
the insulator into a transition to the metallic state), eventually the distance d is exceeded. 
In the limit d, the screening of U{q,u) becomes significant as long as g > l/^'^qxp (see 
Eqs. (0) and (^)). The integration over q in (|1]) is therefore divided into two regions: for 
< g < 1/^'^qTF, U{,q,^) is approximated by the unscreened form, while for q > l/^^g^F, 

egTFsmh(gaj 
The dominant contribution to pt in this regime is 

I h ( kBTV HqTFp , . 

^* ^ 4877 e2 \hDn) {qTrdf ' ^ ' 

Note that in both Eqs. (|^) and (H) pt is an increasing function of ^ (typically g^F ~ 
hence qrp^ ^ !)• However, for ^ ^ d the dependence on ^ is weaker: there is a considerable 
contribution to the frictional drag from density fluctuations on length scales within ^. In 
addition, in this regime the enhancement of screening competes with the enhancement of 
conductivity. 

Finally, when ^ is further increased, eventually the approximations made in the beginning 
of this section break down. For ^ oo, the form of x(g, t^) coincides with that of the diffusive 



state considered in Ref. |TT|. In the low T limit, where <^ h/r, Zheng and MacDonald 
find Pt ~ THn{kBTT/h)/D^. Notice that the T-dependence is altered in the opposite limit 
fc^T ^ h/r: then, in the relevant range of q and u, Imx(g, cj) ~ {dn/ dp)Dq^ / uj, and the 
trans-resistance is 

where i = vpT is the mean free path. When a control parameter is varried so that ^ gradually 
increases from ^ <ti d to ^ ^ oo, pt interpolates between the expressions (|]), (||) and (|). 



III. EFROS SHKLOVSKII INSULATORS 

The effect of Coulomb interactions within the layers on the dissipative processes asso- 
ciated with the hopping mechanism is quite subtle [^. The dominant processes, at small 
but finite uj, q and T, are transitions of electrons between two localized sites that are close 



in energy [^,^]; typically, different pairs of such sites are sparsely distributed. In the pres- 
ence of Coulomb interactions, on one hand the transitions are enhanced due to the greater 
probability to find a singly occupied pair; on the other hand, a Coulomb gap is introduced 
in the single-electron density of states near the Fermi level. The a.c. conductivity a{q,uj), 
assisted by resonant transitions, is more strongly affected by the former, and is therefore 
enhanced compared to the non-interacting case. This has crucial consequences on the drag 
between layers at the ES insulating state. 

Similarly to the previous section, the trans-resistance is evaluated employing Eqs. (P, 
(0) with the appropriate form of assumed to be identical in the two layers. In this 

case the effect of x(g, u) on the inter-layer interaction can be neglected, and U is approx- 
imated by the unscreened form, U{q,uj) ~ V{q)e~'^'^. As will become evident below, the 
tu-integration in Eq. (|I]) has an infra-red divergence, and hence the prominent contribution 
to Pt arises from the lower cutoff. This cutoff is introduced at finite T by dephasing, asso- 
ciated with phonon-mediated relaxational processes, which suppresses resonant transitions 
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between sites separated by a distance r^i larger than the dephasing length L^p. The pair arm 
diverges with the frequency as 



= ^ln{uJo/uj) , uo = (10) 



and is set by the hopping distance 



L,^,m"\ fr..5-l. (u, 



Hence, coherent frequency-driven hopping occur a.t u > Uc where 

u;^ = ^^e'^^l^ = cuoe-^^^/^^'^' . (12) 

To proceed with the calculation of pt using Eq. (|l]), I relate lmx{q,ui) to the a.c. con- 
ductivity through 

lmxiq,uj) = — -a{q,u) . (13) 

a{q, uj) at a finite q was calculated for a 2D ES insulator by Aleiner and Shklovskii who 
obtained 

cr(g, cu) ~ Cie^uj for g < rj^ , 

for rji«g«r\ (14) 

here (j = 1, 2, 3) are numerical constants of order unity, and rj is non-vanishing in case the 
localized single electron states have multi-fractal structure. Note that Ref. ||25| considers 
the limit hu 3> ksT, while in the present case the relevant range of frequencies (which 
dominates Eq. (|l|)) is hu < kBT; however, following Refs. [^] and it can be shown 
that the result differs only by the values of the Cj's. Inserting Eq. ([liD in ([T3| ) yields 
approximate expressions for Imx(g,u;) in three different regimes of q. I then evaluate pt, 
similarly to the previous section, distinguishing the limit cases ^ d and ^ ^ d. 
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In the case ^ -C the high q regime {q > ^) is exponentially suppressed, and the 
g-integration in Eq. (|1])) yields (for carrier density n in the two layers) 

The upper limit = UqC'^'^^^ corresponds to = 2d; the integration over the frequency 
range u > gives a subdominant contribution, that is neglected here. The final expression 
for pt is dominated by the lower cutoff (for sufficiently low T, d): 

clh 1 /r\=^ f/ro\i/2' 



In the opposite regime, where the localization length ^ greatly exceeds d, one should 
account for the contribution of short length-scale density fluctuations with 1/ ^ <^ q <^ 1/ d. 
The integration over q is facilitated by the approximation r^^ ^ d (note that in the relevant 
renge of tu, > ^), and gives 

f^n'e p du f {Cf/Q - CI) {Cl + Ci/2{1-^)) ] 
^* ^ 2n2e2 sinh2(/5;io;/2) \ r^^ r^e j ' ^ 

Similarly to the short ^ limit, the lower cuttoff dominates. The final expression for pt for 
^ > (i is 

{Cl + Cl/2{l-r^)) h 1 /TX^ /r^oV/n 
Z -^-TTI exp — . (18) 



This result essentially differs from Eq. ([T^) only by the algebraic dependence on ^ - short 
length-scale fluctuations are effectively cuttoff by ^ rather than d (the latter sets a higher 
upper cutoff on q). Hence, there is no dependence on the inter-layer separation in this case. 

IV. DISCUSSION AND SUMMARY 

As shown in the calculations detailed above. Coulomb drag between layers in the Ander- 
son insulator state can serve as a clear signature of the presence or absence of a Coulomb gap 
in the layers. Most importantly, the trans-resistance in a double-layer of Mott insulators is 
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suppressed with a decreasing temperature and localization length; in contrast, the presence 
of a Coulomb gap in ES insulating layers leads to a divergence of pt at low T and C,- 

The physical origin of this distinction is the relative enhancement of resonant hopping 
processes in the presence of a Coulomb gap, which ensures a greater probability that a 
"destination site" is empty. In the single-layer a.c. conductivity a{q,uj), this difference is 
indicated by the power-law dependence on u {u versus cu^), though in both cases a vanishes 
at — s> 0. However, the low frequency limit of the dynamical structure factor, 

h 

S{q lmx{q,uj) , (19) 

is crucially distinct: since S ~ cr/'^^, in the Mott insulator S{q,u! 0) is finite, while in 
the ES state it diverges as l/u. As expressed by Eq. (|I|), the frictional drag directly probes 
the dynamics of density fluctuations through the convolution of S{q, u) in the two layers. 
Therefore, in the ES insulator pt depends on the lower frequency cutoff, associated with the 
dephasing length, hence diverges at low T and ^. In the Mott insulator this anomaly does 
not exist, and pt decreases with ^ due to the suppression of excited density fluctuations - in 
agreement with the naive intuitive argument. 

In both types of insulators, pt{0 depends on whether ^ is smaller or larger than the 
layers separation d. The predictions of this paper can be summarized as follows: (a) in the 
Mott insulator, pt ~ similarly to the free electron gas case. The localization is manifested 
as a strong dependence on ^ for ^ ^ d, 

Pt oc {^/df , (20) 

while foTC,^d the subtle role played by screening effects yields a weaker dependence, 

Pt (X HqTpO ■ (21) 
(b) in the ES insulator, the resulting pt{C,,T) is 

''^^^©''''^{©'^l ^^^^ 
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where ksTo = e^/e^. Note that pt diverges at low T and ^ with the same exponential factor 
as the in-layer resistance p. The algebraic prefacor typically suppresses pt with respect to 
p by 3 to 4 orders of magnitude. This ensures that although the drag is a huge effect in 
this case, the two-layer resistivity tensor is still almost diagonal, and the weak-coupling 
assumption underlying Eq. (|I|) is justified. 

Experimentally, a crossover from an ES to a Mott insulator can be in principle controlled 
by a metallic back gate, which effectively attenuates the range of interactions in the layer. 
The above predictions imply that at the ES state, the onset of an insulating behavior should 
be accompanied by a sharp increase of pt towards the insulating regime; when the insulating 
state is well approximated by the Mott behavior, pt should be peaked near the transition 
to the insulator, and strongly attenuated when ^ is reduced below d. A similar qualitative 
behavior is also expected in case where some of the simplifying assumptions of this paper 
are violated - for example, when the layers are not identical, and in particular if only one 
of them undergoes a transition to the insulator. However, note that in the case where one 
of the layers is a good conductor, it may serve as a back gate which suppresses Coulomb 
interactions in the other. 

In order to observe an appreciable drag at low T in the Mott state, the desired experi- 
mental setup should involve low mobility, low density samples. For £ ~ 1 pm, n ~ 10^° cm~'^, 
d = 200 A, band mass of m ~ 0.07me, e = 13 (typical to GaAs) and T = 1 i^, I get pt in the 
order of a few tens of rriQ in the regime > d. When the localization length is reduced to 
20 A, Pt is attenuated by a factor of 10^^. Assuming the same parameters in the ES state, I 
estimate pt ~ 575 for ^ ~ 1000 A; at ^ ~ 100 A, the trans-resistance rises to pt ~ 100 kO,. 

Finally, the experimental testing of the effects predicted in this paper involves a number 
of difficulties. Primarily, once the in-layer resistance is comparable to that of the barrier 
separating the two layers, tunneling across the barrier is no longer negligible; its contribution 
should be carefully eliminated. In addition, to obtain a sizable signal, the voltage imposed 
on the drive layer in an insulating state may be large enough to produce non-linear response 
effects. The role of inter-layer coupling mediated by phonons (see, e.g., Gramila et al. in 
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Ref. 0, and Ref. [0), and possible thermoelectric effects (Solomon et al. in 0, Laikhtman 
et al. in [0), are not discussed in the present paper. These are generally expected to be 
subdominant at low T; however, at the onset of an insulating state (and particularly near a 
MIT) such effects may be enhanced as well. Nevertheless, the estimates made above indicate 
that Coulomb drag near the onset of an insulator is an appreciable effect, and a sensitive 
probe of the significance of Coulomb interactions in insulating states. 
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